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GEOMETRIC  IJiTERPRETAT  I  ON  OF  LANTZIG  S  CONVEX  PROGRAMMING  ALGORITHM 

by  Pierre  Huard 


I.  Review  of  the  Algorithm  [1] 

We  consider  the  following  program: 


(1) 


PI 


jTind  x  e  R  such  that  0Q(x)  is  maximum 
'-subject  to  01(x)  ^  0  for  i  =  1,2, . . .  ,m  . 


where  0Q  ,  01  are  concave  in  x  and  R  is  compact  and  convex. 

We  replace  PI  by  a  linear  program  PII-K  which  gives  a  lower 
approximate  solution  0q(x)  for  PI.  This  linear  program  is  constructed 
from  a  finite  number  of  points  x  =  X^  ,  j  e  K,  which  belong  to  R  in 
the  following  way: 


^ind  Xj  0  ,  j  e  K  such  that 


I 


j€K 

PII-K  /  subject  to 


0Q(X^)Xj  is  maximum 


(2) 


(3) 


JSK 


Iv 


j€K 


1,2, * « • , m 


l£t  0  be  a  column  vector  whose  components  are  0  ,  i  =  1,2, _ ,m  . 

=  (0(X^)  |  j  e  K)  ,  matrix  of  columns  0(X^) 

f  =  (0O(XJ)|  j  £  K)  ,  row  vector  of  components  0Q(X^) 

K  1 

e  =  (1,1, ...,l),  row  vector  of  components  e^  =  1  ,  j  e  K 
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The  linear  program  PI1-K  and  the  corresponding  Kuhn-Tucker  conditions 
in  matrix  notation  read: 


Find  0  and  maximize  fK^K 

PII-K  subject  to 
(2)  £  0 
(3)  e1^  =  1 

i 

\ 

w 

(5) 

Kuhn- 

(6 )  Tucker 
conditions 

(7) 

The  dual  variables  Jl(vector)  and  n°(scalar)  correspond  respectively 
to  the  relations  (2)  and  (3). 

The  optimization  of  the  linear  program  PII-K  yields  an  optimal 

solution  ,  at  finite  distance,  and  corresponding  dual  variables  fi 

-0 

and  H  as  functions  of  K  . 

One  solves  then  the  following  auxiliary  program 

Find  x  €  R  such  that 

PIII-K 

m 

^ni0i(x)  +  0q(x)  is  maximum 
i=l 

let  x  =  Xs  be  the  optimal  solution  of  PIII-K  .  The  index  s  is  then  added 
to  the  set  K  of  the  indices  which  define  the  linear  program  PII-K  . 


n  (n°)  ^  0 

hak  +n°eK  +  fK  1  0 

ha^  =  o 

[ITaK  +  n°eK  +  fK]  \  =  O 


PII-K  has  then  one  more  variable  A  .  If  this  new  variable  is  a 

s 

candidate  for  entering  the  basis,  then  PII-(K+s)  can  be  improved, 
if  not  then  PII-K  yields  the  optimal  solution  x  for  PI  ,  viz: 

(8) 

j£K 

or  with  obvious  notations 

(9)  x  =  x\  • 

If  the  algorithm  is  infinite,  and  with  the  restriction  of  a  nondegenerate 
solution  to  the  initial  PII,  the  solutions  of  PII-K  converge  to  the 
optimal  solution  of  PI. 

II.  Geometric  Interpretation  of  PII 

If  we  consider  the  point  defined  by 


(10) 

x(V  = 

where  y? 

has 

the  same  meaning  that  in  (9),  we  have  for  ell  feasible 

solutions 

of  PII: 

(11) 

0Cx(XK)]  ^  a\.  £  0 

(12) 

0oMy^  f\ . 

In  fact,  AK  A^.  and  f^A^  rePresen'fc  barycentric  interpolations 
of  concave  functions.  From  relation  (ll)  we  know  that  to  every  feasible 
solution  A^.  of  PII-K  corresponds  a  feasible  solution  X(^K)  bo  PI, 
i.e.,  that  the  domain  of  the  solution  x(A^)  obtained  from  feasible  solu¬ 
tions  A^.  for  PII  is  contained  in  the  domain  of  feasible  solutions  of 
PI.  Moreover,  relation  (12)  shows  that  the  value  of  the  objective  function 


- k - 


III.  Geometric  Interpretation  of  PHI 
PHI  has  the  following  form: 

Find  x  e  R  such  that  It0(x)  +  0  (x)  is  maximum 

PHI 

.where  R  £>  o  is  a  given  vector. 


We  have  seen  that  PHI  is  used  to  determine  the  candidate  variable  of  the 
simplex  method  to  enter  the  basis  of  PII  (linear  program  in  A). 

In  the  standard  simplex  method,  each  step  corresponds  to  an  extreme 


point.  To  each  edge  leaving  that  vertex  corresponds  a  non -basic  variable 

A.  .  We  then  determine  an  edge  leaving  that  vertex  which  gives  the  greatest 
J 

increasing  slope  of  the  objective  function.  The  slope  being  taken  with 


respect  to  the  nonbasic  variable  A.  corresponding  to  that  slope.  We  can 

J 

also  select  any  edge  yielding  a  strictly  positive  slope  (not  necessarily 
the  largest). 


Here  the  problem  is  different:  we  do  not  have  these  edges,  but  we 
have  to  construct  one,  with  positive  slope  (or  the  largest  positive  slope), 
increasing  the  A-space  by  one  dimension.  The  position  of  the  new  edge  and 


the  position  of  the  new  gradient  to  the  objective  function  are  related  and 
are  functions  of  the  selection  of  a  point  x  =  Xs  of  R  .  On  the  contrary, 


the  projections  of  the  normals  on  the  new  faces  of  PII  as  well  as  that 
of  the  new  gradient  In  the  old  A-space,  remain  constant. 
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Fig-  3 

When  PIII-K  has  generated  a  candidate  ,  the  first  iteration  of 
(K+s )  is  the  introduction  of  in  the  basis,  because  it  is  the  only 

candidate.  But,  once  the  basis  is  changed,  there  may  exist  other  candidates, 
namely  slack  variables.  Note  that  PIII-K  yields  an  additional  variable  and 
candidate  for  entering  the  basis  to  FII-K  ,  the  criterion  of  choice  of 
that  candidate  takes  only  into  account  the  relative  gain  (or  slope,  as  we 
said  above)  of  the  first  iteration  of  PII-(K+s)  and  not  the  total  gain 
of  the  following  iterations  necessary  for  the  optimization  of  PII-(K+s)  .  ' 
The  optimization  of  PIII-K  does  not  necessarily  yield  the  "best" 
candidate  for  PII-(K+s),  in  the  sense  of  improvement  for  the  objective 
function  of  PII-(K+s).  Let  us  illustrate  this  by  an  elementary  example. 
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EXAMP LE  Let  x  be  a  scalar,  M  a  very  large  scalar. 


con- 


We  have  also  d  =  (classical  result  of  the  parametrizatlon  of  the 

stant  column  of  a  linear  program) 


(14) 

(15) 


\  "s 


_i  = 

al  \  +  \ 


=  * 


(similar  triangles) 

(similar  triangles)  . 


It  follows  that  the  value  of  the  slope  Is: 

(16)  p"f'  +  TC“  nS.  +  *o  ' 

s  s 


This  slope  Is  a  function,  "by  a.^  and  aQ  ,  of  x  =  Xs  given  by  FIII-K  . 
Letting  x  be  the  abscissa  of  the  point  (J) : 


(17)  ax(x)  =  01(x)  -  -i  (x  -  x)  +  0 

(18)  aQ(x)  =  0Q(X)  +  f  (x  _  x)  +  C 


tt  follows 

(19)  p(x)  =-n1a1(x)  +  aQ(x)  =  n101(x)  +  0Q(|  -  ^)(x  -  x) 

=  It  0x(x)  +  0Q(x)  +  C 

We  see  that  for  the  objective  function  of  PHI  the  part  ir*"0^(x)  has 

for  origin  the  term  d/X  ,  i.e.,  the  increase  of  the  optimal  value  of 

s 

the  objective  function  of  PII  due  to  the  enlargement  of  the  domain  of 
feasible  solutions,  whereas  the  part  0g(x)  has  the  term  b0/>>s  for 
origin,  corresponding  to  a  better  approximation  of  the  function  0Q(x) 
of  PI  for  the  new  form  of  the  objective  function  of  PII. 

These  results  can  be  reproduced  for  the  n-dimensional  case. 
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Let: 


be  an  optimal  solution  of  PII-K 


We  have : 


XjJ.  the  solution  obtained  after  the  first 
iteration  of  PII-(K+s) 

I  the  optimal  basis  of  PII-K,  with  respect 
bo 


E 


the  set  of 
of  PII-K 
for  X  =  ^ 


indices  of  the  constraints  (2) 
satisfied  exactly  (equality) 


B  = 


basic  matrix  of  PII-K 


(xK,  X6) 


Variation  in  the  objective  function  of  PII  from 
x  to  x' 

=  Variation  of  the  old  objective  function  of  PII 
from  x  to  x' 

+  Variation  at  the  point  x 1  when  we  replace  the 
old  objective  function  by  the  new  one. 


More  precisely: 


Variation  of  the  old  function  from  x  to  x'  = 

(20) 

=  f1^’  -  xI) 


=  (HAS  +  n°)X' 
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THE  PROBLEM  OF  THE  MAXIMUM  FLOW  IN  TRANSPORTATION 
WITH  CORRESPONDING  CONSTRAINTS 

A.  Ghouila-Houri 

(Member  of  the  Operations  Research  Group  of  the  S.  N.  C.  F.,  France) 

In  a  transportation  network  in  which  for  each  entry  arc  there  corres¬ 
ponds  an  exit  arc  and  vice  versa,  corresponding  pairs  of  entries  and  exits 
are  required  to  carry  the  same  flow.  Under  these  conditions,  we  can  obtain 
the  maximum  flow  by  applying  the  Ford-Fulkerson  algorithm  to  a  sequence  of 
networks,  each  of  which  is  obtained  from  the  labeling  of  the  preceding  net¬ 
work.  The  final  result  is  attained  in  a  finite  number  of  steps. 

1.  Introduction 

The  study  of  the  problem  of  engine  scheduling  at  the  S.  N.  C.  F.  led 
us  to  the  following  theoretical  problem: 

We  consider  a  transportation  network  in  which  there  is  a  two-way 
correspondence  between  entry  and  exit  arcs. 

COLLECTION  OF  NODES 
OTHER  THAN  SOURCE 
AND  SINK 


FIG.  I 
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To  each  arc  u  in  the  network  is  attached  a  capacity  C(u).  The 
problem  is  to  send  flow  in  this  network  subject  to  the  capacity  con¬ 
straints:  0(u)  £  C(u ) ;  the  correspondence  constraints:  flow  through  an 
entry  arc  =  flow  through  the  corresponding  exit  arc,  and  such  that  the 
flow  is  maximum. 

2.  The  Problem  of  the  Maximum  Flow  in  a  Transportation  Network  (To 
review  known  notions) 

We  shall  recall  the  principles  of  the  Ford-Fulkerson  algorithm 
which  allows  us  to  solve  the  following  simpler  problem: 

Given  a  transportation  network;  for  example,  the  one  shown  in 
Fig.  2,  in  which  the  numbers  between  parentheses  represent  arc  capa¬ 
cities,  the  problem  is  how  to  send  maximum  flow  from  source  to  sink. 


A  t;t  B 


FI  G  2 


The  algorithm  of  Ford-Fulkerson  consists  of  the  following  operations 
(l)  Send  an  arbitrary  flow:  as  illustrated  in  Figure  3. 
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FI  G.  3 


(2)  Proceed  with  labeling  nodes  according  to  the  following  rules: 

a)  S  must  be  labeled  (with  X) 

b)  For  a  non-saturated  arc  (x  ,  y)  if  x  is  labeled,  label  y 

c)  For  a  non-empty  arc  (y  ,  x)  if  x  is  labeled,  label  y 

(3)  a)  if  S'  is  labeled,  we  can  improve  the  flow  by  one  unit.  For 

example,  in  Figure  5,  S  is  labeled  by  rule  a;  C  ,  F  ,  and  G 
by  rule  b;  E  by  rule  c;  B  by  rule  b;  A  by  rule  cj  D  and 
S  by  rule  b.  If  we  add  one  unit  of  flow  to  SC,  CF,  FG,  sub¬ 
tract  one  unit  of  the  flow  from  EG,  add  one  unit  to  EB,  subtract 
one  from  AB,  add  one  unit  to  AD  and  DS'  we  obtain,  as  shown 
in  Figure  4,  a  one  unit  improvement  over  the  previous  flow. 
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FI  G.  4 

b)  If  it  is  not  possible  to  label  S'  (for  example  in  Figure  4, 
only  S  can  be  labeled),  then  maximum  flow  is  attained.  To 
prove  this,  consider  the  set  of  arcs  U+  such  that  the  be¬ 
ginning  nodes  are  labeled  and  the  end  nodes  are  unlabeled, 
and  the  set  of  arcs  U~  such  that  the  end  nodes  are  unlabeled 
and  the  beginning  nodes  are  labeled.  The  flow  through  the  net¬ 
work  from  S  to  S 1  is  also  the  flow  out  of  the  set  of  labeled 
nodes  (containing  S)  to  the  set  of  unlabeled  nodes  (containing 
S').  Accordingly,  it  is  the  sum  of  flow  in  arcs  going  from  the 
first  set  to  the  second  minus  the  sum  of  flow  going  from  the 
second  to  the  first  set 
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i  =  l,  2,  ...  m  and  where  -»  <  p  <  oo  .  The  arcs  of  this  graph  and 
their  capacities  will  be  defined  in  the  following  manner: 


I  31  ENTRY  ARC 


nd  /  <2> 

2  Entry  arcx 


X1  i,> 

V  , 

S.  /  *  /  1X1 


iST  exit  arc 


2ndEXIT  ARC 


INITIAL  NETWORK 


x'p 


4'XaSSC\mJBtJt> 


-12>  -/  (i)  d/ 


— 


l'/*p  ''I  /r  \ — '= 1  ^ 


AUXILIARY  GRAPH 


FI  6.  5 


a )  For  each  value  of  p  ,  the  nodes  x1  ,  x2  , .  . . ,  x"1  are  con 

P  P  P 

nected  by  arcs  exactly  like  the  one  in  the  initial  network. 

b)  For  each  entry  or  exit  arc  in  the  initial  network  say  Sx1 

x^S'  ,  and  for  each  p  there  corresponds  an  arc  x^x*  ,  w: 


a  capacity  equal  to  the  smallest  of  the  two  capacities. 

Figure  5  illustrates  the  construction  of  such  an  auxiliary  graph. 

In  a  network,  send  flow  which  assigns  to  every  arc  u  a  positive 
integer  flow  n(u)  £  C(u)  such  that  at  every  node  other  than  source 
and  sink  the  conservation  equations  are  satisfied;  i.e. ,  the  sum  of 
the  flows  in  entering  arcs  equals  the  sum  of  the  flows  in  exit  arcs. 

We  say  likewise  for  the  auxiliary  graph  G  that  we  have  sent  a  flow  into 
G  if  to  each  arc  u  is  assigned  a  positive  integer  flow  n(u)  C(u)  so 
that  at  each  node  the  conservation  equations  are  satisfied.  We  say  that 
the  flow  is  periodic  if  the  flow  is  the  same  in  any  two  homologous  arcs, 
where  homologous  arcs  are  two  arcs  which  correspond  to  the  same  arc  or 
pair  of  arcs  in  the  initial  network.  Under  these  conditions,  it  turns 
out  to  he  the  same  to  send  flow  in  the  initial  network  obeying  the  corres¬ 
pondence  constraints  or  sending  periodic  flow  in  the  auxiliary  graph  G. 

U.  Expression  for  the  Quantity  to  he  Maximized 

The  quantity  to  be  maximized,  which  in  the  initial  network  is  the 
sum  of  flows  in  entering  arcs,  is  found  in  the  graph  G  as  the  sum  of 
flows  in  the  arcs  connecting  nodes  with  subscript  p  to  the  nodes  with 
subscript  p  +  1  .  However,  we  can  give  it  a  more  general  expression. 
Consider  a  function  p(i)  which  assigns  to  every  value  of  the  super¬ 
script  i  a  value  of  the  subscript  p  .  Each  function  p(i)  defines  a 
partition  of  nodes  of  graph  G  into  two  sets:  the  set  E[p(i)]  of  x* 
such  that  p  p(i)  ,  and  the  complementary  set.  Calling  U[p(i)]  the 
set  of  arcs  for  which  only  the  origin  node  belongs  to  E[p(i)]  ,  W[p(i)] 
the  set  of  arcs  for  which  only  the  terminal  node  belongs  to  E[p(i)]  and 
finally,  given  a  periodic  flow  passing  through  G  ,  we  call  n(u)  the 
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flew  in  arc  u  .  If  we  take  for  the  p(i)  the  constant  function  of 
i  ,  say  p(i)  =  pQ  ,  then  the  quantity  to  he  maximized  according  to 

what  we  mentioned  previously  is 


£  n(u)  ; 
ueU[p(i)) 

on  the  other  hand  for  such  a  function  constant  in  i  ,  the  set  W[p(i)] 
is  empty  since  we  did  not  have  in  the  construction  of  G  arcs  with  the 
origin  node  having  subscript  p  +  1  ,  and  terminal  node  having  subscript 
p  .  Accordingly,  we  can  put  the  quantity  to  be  maximized  in  the  form 

z 

ueu[p(i)]  uew[p(i)] 

where  p(i)  is  a  function  constant  in  i  .  From  the  flow  conservation 
equations,  we  can  easily  deduce  that  the  quantity 


n(u)  - 


z 


n(u) 


^  n(u)  -  n(u) 

u€U[p(i)]  u£M[p(i)] 

does  not  depend  on  p(i)  and  we  denote  this  quantity  by  F[n'u)]  . 

If  we  define  C[p(i)3  =  £  C(u)  the  following  theorem  is  obtained 

ueu[p(i)l 

THEOREM:  For  any  function  p(i)  ,  F[n(u)]  ^  C[p(i)]  .  Indeed  0  £  n(u) 

£  C(u)  and  hence 

F[n(u)]  =  V  n(u)  -  £  n(u)  £  £  C(u)  =  C[p(i)]  . 
u€U[p(i)]  u€W[p(i)]  ueu[p(i)3 
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The  nodes  of  G[p(i)]  are  the  elements  of  A[p(i)]  = 

E[p(i)]  -  E[p(i)  -  1]  to  which  are  added  four  points  I  ,  J  , 

I  ,  J  . 
o'  o 

The  arcs  of  G[p(i)]  are  the  arcs  of  G  after  the  following 
transformation:  (i)  replace  by  I  any  node  of  E[p(i)  -  1]; 

(ii)  replace  by  J  any  node  not  in  E[p(i)];  (iii)  discard  any 
arc  with  both  ends  at  the  same  node;  (iv)  add  the  arcs  (i  ,  I) 
and  (j  ,  Jq)  both  with  infinite  capacity. 

Three  categories  of  arcs  will  be  distinguished: 


The  set  of  arcs  which  do  not  touch  I  or  J  ; 

the  set  Ug  of  arcs  whose  origin  is  I  or  whose  end  is  J  ; 

the  set  of  arcs  whose  end  is  I  or  whose  origin  is  J  . 

Figure  6  shows  a  network  G[p(i)]  corresponding  to  the  example  of 
Figure  5  taking  p(l)  =  p(2)  =  p(5)  =  p(4)  =  1,  p(5)  =  p(6)  =  p(7)  =  0  . 


I«t  N[p(i)]  be  the  value  of  the  maximum  flow  in  G[p(i)].  Then 
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N[p(i)]  is  nonnegative  since  a  zero  flow  is  always  possible.  In 

addition,  N[p(i)]  £  G[p(i)]  where  C[p(i)]  is  the  capacity  of  the 

arc  consisting  of  the  arcs  other  than  J  ,  having  one  end  at  J  . 

The  Ford-Fulkerson  algorithm  applied  to  this  network  yields  a 

partition  of  the  nodes  of  G[p(i)]  into  labeled  nodes  (including  I  ) 

o 

and  unlabeled  nodes  (including  Jq)  .  For  p(i)  let  0p(±)(i)  be  a 

corresponding  function  equal  to  p(i)  if  x1  is  labeled,  to 

p(i)  -  1  if  is  not  labeled.  Then  C  C0p  ^ ±  j C 1 ) ]  is  precisely 

the  capacity  of  the  cut  indicated  by  the  preceding  partition,  so  then 
N[p(i)]  =  C[0p(i)(i)]  . 

7.  Solution  of  the  Given  Problem 

Suppose  for  the  moment  that  p(i)  is  such  that  N[p(i)]  = 

C[p(i)]  and  the  flow  in  G[p(i)]  is  maximum.  Then  the  arcs  of  U2 
are  saturated,  and  those  of  are  empty. 

Furthermore,  an  arc  of  and  an  arc  of  U  for  A  ^  p  or  for 

A  =  p.  =  1  cannot  correspond  to  identical  or  homologous  arcs  of  G 
and  any  arc  u  of  G  is  homologous  to  at  least  one  arc  corresponding 
to  an  arc  v  of  G(p(i)]  .  If  u  e  ^  ,  let  n(u)  be  the  flow  in  v  . 

If  v  e  U2  ,  let  n(u)  the  flow  in  v  .  If  v  e  U,  ,  let  n(u)  =  0  . 

5  J 

A  periodic  flow  in  G  is  obtained  in  this  way.  The  periodicity 
is  trivial,  and  it  is  easily  seen  that  the  capacity  constraints  and 
conservation  equations  are  satisfied.  This  flow  is  maximum  because 

F[n(u)]  =  ^  n(u)  -  ^  n(u)  =  ^  C(u) 

ueU[p(i)]  u€W[p(i)]  ueU[p(i)] 

-  \  0  =  C[p(i)] . 

u€W[p(i)] 
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and  a  previous  theorem  shows  that  any  flow  in  G  satisfies  F[n(u)]  £ 
C.Ip(i)]  . 

Hence,  we  can  obtain  the  answer  to  the  given  problem  once  a  func¬ 
tion  p(i)  is  known  such  that  N[p(i)]  =  C[p(i)]  . 

Thus,  the  following  theorem  is  important: 

THEOREM:  There  exists  p(i)  such  that  N[p(i)]  =  C[p(i)]  . 

PROOF:  Consider  the  sequence  of  functions  pQ(i)  ,  p.(i), . . .,Pn(t), . . . . 

defined  beginning  with  the  arbitrary  function  PQ(i)  by  the  recurrence 

relation  Pn+j_(i)  =  ^  \  ( i )  •  To  it  corresponds  a  sequence  of  posl- 

Pnv  / 

tive  integers  C[po(i)]  ,  C[p1(i)], . . . ,  C[pn(i)],...  .  This  sequence 
cannot  be  strictly  decreasing.  Hence,  there  is  an  integer  n  such  that 

CtPn(i»  -  CtPn+iC1)]  •  But  N[pn(i)]  =  C[0p  (±)(i)]  =  C[pn+1(i)]  . 

Hence  N[pn(i)]  =  C[pn(i)]  . 

The  preceding  proof  defines  an  algorithm  which  will  be  applied  to 
the  network  of  Figure  ^  as  an  example: 

To  begin,  let  us  take  PQ(i)  =  0  .  G[pQ(i)]  is  the  initial  network 
(except  for  the  modifications  of  the  capacities  of  entry  and  exit  arcs). 
Figure  7(a)  shows  a  maximum  flow  in  this  network  (with  the  same  flow  in 
corresponding  entry  and  exit  arcs). 


FIG.  7a 
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A  REVIEW 


A  TECHNIQUE  FOR  RESOLVING  DEGENERACY  IN  LINEAR  PROGRAMMING 

by 

Philip  Wolfe 

J.  Soc.  Indust.  Appl.  Math.  11  (1963),  205-211 

Prepared  by  George  B.  Dantzig 
for 

"Mathematical  Reviews" 

Geometrically,  the  simplex  method  for  solving  linear  programs 
passes  iteratively  from  one  extreme  point  to  a  selected  neighbor  in 
a  convex  polyhedral  set.  If  the  set  is  defined  by  a  system  of  m 
linear  equations  in  n  non-negative  variables  then  algebraically 
an  extreme  point  solution  E  is  obtained  if,  setting  n  -  m 
variables  equal  to  zero,  there  is  a  unique  non-negative  solution  in 
the  remaining  "basic"  set  of  m-variables.  If  the  same  extreme  point 
E  corresponds  to  more  than  one  basic  set,  the  solution  is  called 
"degenerate."  If  so,  the  value  of  linear  "objective"  form  (to  be 
minimized)  may  not  decrease  from  iteration  to  iteration.  Several 
special  selection  rules  for  successively  choosing  basic  sets  have 
been  proposed  to  guarantee,  under  degeneracy,  termination  in  finite 
number  of  steps. 

According  to  the  author,  his  paper  "is  closely  related  to  the 
material  of  Dantzig 's  Inductive  proof  of  the  simplex  method.  That 
paper  may  be  viewed  as  showing  the  existence  of  a  class  of  choice 
rules  which  prevent  cycling  in  the  simplex  method  and  the  present 
paper  viewed  as  exhibiting  a  member  of  that  class."  At  any  iteration 
let  R  be  the  set  of  equations  in  the  simplex  tableau  whose  constant 
terms  are  zero  and  from  which  the  pivot  term  can  be  selected.  These 


constant  terms  are  each  replaced  by  the  polynomial  0  +  e  where  e  >  0  . 
(In  selecting  pivots,  polynomial  expressions  in  e  and  later  e2,... 
are  compared  using  a  lexicographic  ordering  of  their  leading  coeffi¬ 
cients).  Degeneracy  becomes  less  "deep''  if  a  pivot  term  occurs  outside 
of  R  ;  if  so,  there  will  be  a  positive  decrease  of  the  objective  form. 
However,  degeneracy  deepens  if  a  proper  subset  of  R'  of  the  R  equa¬ 
tions  should  develop  polynomial  expressions  in  e  with  all  zero  coeffi¬ 
cients.  For  these  equations,  the  constants  are  replaced  by  0  +  Oe  +  e2 
and  pivots  are  now  selected  from  R1  ,  etc. 

The  author's  proof  is  essentially  inductive  (although  not  exactly 
so  stated).  The  number  of  pivots  selected  in  the  subset  R  will  be 
finite  because,  by  induction,  there  can  be  only  a  finite  number  of  steps 
in  R'  which  results  in  an  optimum  solution  or  a  non-zero  decrease  in 
the  €  term  of  the  objective  form.  This  implies  nonrepetition  in  the 
choice  of  basic  sets  of  variables. 

The  paper  concludes  with  seme  empirical  observations  on  the  number 
of  iterations  required  to  solve  some  simple  degenerate  problems  using, 
unfortunately,  a  different  pivot  selection  rule. 
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